Origin and Transport Signatures of Spin-Orbit Interactions in One- and 
Two-Dimensional SrTiOs-Based Heterostructures 
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We study origin of Rashba spin-orbit interaction at SrTiOs surfaces and LaAlOa/SrTiOa interfaces 
by considering the interplay between atomic spin-orbit coupling and inversion asymmetry at the 
surface or interface. We show that, in a simple tight-binding model involving 3d t2g bands of Ti 
ions, the induced spin-orbit coupling in the d^z and dy^ bands is cubic in momentum whereas the 
spin-orbit interaction in the dxy band has linear momentum dependence. We also find that the spin- 
orbit interaction in one-dimensional channels at LaAlOs/SrTiOa interfaces is linear in momentum 
for all bands. We discuss implications of our results for transport experiments on SrTiOs surfaces 
and LaAlOs/SrTiOa interfaces. In particular, we analyze the effect of a given spin-orbit interaction 
term on magnetotransport of LaAlOg/SrTiOs by calculating weak anti-localization corrections to 
the conductance. 



I. INTRODUCTION 



The metallic interfaceaS between the insulators 
LaAlOa (LAO) and SrTiOa (STO) exhibit both 
superconductivity^! and magnetisrrPEl a,nd, for some car- 
rier densities, they occur simultaneously^lt^. Moreover, 
there is evidence that there is significant spin-orbit inter- 
action (SOI), as welpHSl. If all three of these phenomena 
are present, then the ingredients are in place for topolog- 
ical superconductivity that could support Majorana zero 
modes in confined structureanHHl. Furthermore, strong 
SOI could pave the way towards spintronicsiSl applica- 
tions of devices based on oxide interfaces. In this paper, 
we give a simple microscopic understanding of the SOI 
at LAO/STO interfaces induced by the combined effects 
of atomic spin-orbit coupling and the interfacial electric 
field. 

Caviglia et alW^ and Ben Shalom et al^ made im- 
portant experimental progress towards understanding 
SOI effects at LAO/STO interfaces. They found evi- 
dence that the magnetoconductance of LAO/STO inter- 
faces could be interpreted as resulting from weak anti- 
localization (WAL). The SOI that they deduced showed 
strong dependence on gate voltage, peaking at or near 
the gate voltage at which the superconducting T^ is max- 
imized. Nakamura et ali^ measured magnetoconduc- 
tance at the surface of STO and found that it could be 
fitted to a cubic Rashba SOI. Zhong et a/.E3 performed a 
density functional theory (DFT) calculation, from which 
they derived an effective tight-binding Hamiltonian. A 
key ingredient supplied by the DFT calculation is the 
magnitude of inter-orbital hopping terms. They used the 
resulting tight-binding Hamiltonian to deduce a Rashba- 
type energy splitting between the two spin components 
of « 2 meV in the d^y band, and a much larger split- 
ting sa 20 meV at the crossing point of the dxy and d^z 
bands. Khalsa et al}^ further elucidated this by showing 
that inter-orbital hopping terms are due primarily to the 
polar lattice displacement at the interface. 



In this paper, we give a simple analysis of the Rashba 
SOI in an effective tight-binding Hamiltonian for the t2g 
bands of STO surfaces and LAO/STO interfaces, in or- 
der to better understand its basic qualitative features. 
Within our effective model, we find that Rashba SOI is 
linear in momentum in the dxy band, in agreement with 
Refs. [TT] and 1181 ^^nd is cubic in momentum in dxz and 
dyz bands. 

We also consider one-dimensional(lD) channels at the 
LAO/STO interface. Cen et al}^ have fabricated such 
channels by 'drawing' them with an atomic force micro- 
scope (AF M) tip . Superconductivity is observed in these 
ID channela^SEi]. Fidkowski et al^ proposed a theory 
for magnetism and superconductivity in such channels in 
which conduction electrons in the channels interact with 
localized spins to catalyze magnetic order and interact 
with local superconducting fluctuations in STO to stabi- 
lize quasi-long-ranged superconducting order. (For other 
theoretical perspectives, see Ref. [^5H^ ) In the presence 
of strong SOI, such superconductivity can support Ma- 
jorana zero modes at the ends of wires. Conduction in 
such a channel will be dominated by mobile dxz or dyz 
electrons, depending on the direction of the channel. We 
show the SOI will be of linear-in-fc Rashba form in this 
case due to the broken rotational symmetry in the plane. 

Based on the aforementioned conclusions regarding the 
different forms of the Rashba SOI in various geometries, 
we compute WAL correction to the magnetoconductance 
and show how one can distinguish different forms of 
Rashba SOI in transport experiments. We discuss the 
relevance of these results to understanding experiments 
at STO surfaces and LAO/STO interfaces. In the latter 
case, the WAL signal can distinguish between transport 
dominated by the dxy band; the dxz,yz bands; or ID chan- 
nels. The ID case should apply to channels 'drawn' with 
an AFM tip as well as to the ID channels that appear 
to occur in putatively two-dimensional (2D) systems^^ al- 
though, in the latter case, it will also be important to 
account for the coupling between different ID channels. 

The paper is organized as follows. In Sec. ITT) we intro- 
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FIG. 1. Dispersion of Ho for Ab = 50meV 



duce our effective tight binding model. The calculation 
of the effective spin-orbit interaction based on this model 
is presented in Sec. |III| In Sees. |IV| and [V] we calculate 
weak anti-localization corrections to the conductance and 
discuss the manifestations of SOI in the context of recent 



experimenti 



II. THREE-BAND MODEL 



The Fermi energy at STO surfaces or STO-based 



interfaces lies in the 3d t2g bands of Ti ions near 
the surface/interface. The t2g bands at the sur- 
face/interface are confined in z-direction, which is normal 
to the surface/interface; consequently, they form a two- 
dimensional electron gas (2DEG). Here we consider only 
the top layer of STO and its t2 a ban ds. The Hamiltonian 
for these bands takes the forrrp32Sl 



H^Ho + Ha 



so 



H„ 



(1) 



Here, Hq is comprised of nearest-neighbor hopping and 
on-site interaction terms that are diagonal in orbital 
space. In momentum space, it can be written in the 
form Hq = Hq (E) <j'^, where a'^ is the identity matrix in 
spin space, and 
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Here A^; is the energy difference between the dxy band 
and the dxz , dyz orbitals due to the confinement along 
z-direction. A recent DFT calculation^ suggests mi — 
0.41?7ie and m/j = 6.8171,, for bulk STO. The second term 
in Eq.([I]) corresponds to the atomic spin-orbit coupling 



term which can be written in the form: 



Haso = 



Zg^Be 



Zfln L ■ a 



L-a 



Vaso- 



(3) 



where the dimcnsionful prefactor Vaso = ^^% '^ — jj 

L = r xp and a = 2S /h with Z being the effective nuclear 
charge on the Ti atoms. The effective nuclear charge for 
the d-orbital electron in a neutral Ti atom Z w 8.1. 

Atomic spin-orbit coupling projected to i2g orbital 
bands can be treated as an on-site orbital mixing term. 
Indeed, let's consider the limit mi,mh ^ od and com- 
pute matrix elements of the Hamiltonian Haso between 
different orbital states: 

(j, dxz, cr'\HASo\j, dxy, cr) = 

Vaso— ^{j, dxz, cr'\-^\j, dxy, cr) (4) 

where \j, dxy, cr) represents a state of an electron of spin 



a in the dxy orbital on site Vj. Given that d. 



<xz ancL dxy 



orbital wavefunctions are both odd in x, the matrix el- 
ements vanish by symmetry: {j,dxz,<j'\Ly\j,dxy,(j) = 
{j, dxz, '^'\Lz\j, dxy, cr) = 0. The non-zero matrix element 
involves dxz and dxy bands 

{j, dxz, >^'\Haso\J, dxy. A) 

^VASO^^{j,dxz,y\ —^\j,dxy,\) 



= VASo[(^x\\',\^-{j,dx 



ypz 



ZPy 



\j, dxy 



^iAAso [crx]x',\- (5) 

In the last line, we have introduced the energy Aaso'- 

AAso = VAsofiZ) (6) 

where the dimensionless form factor f{Z) is defined as 



')^# 


0; 


dxz\ 


-zpy 


\j. 


2Z^ 1 


00 
-00 


dx dy dz 


xze" 
r3 


Zt 
3 


81V / 




Z4 










405' 











dz dy 



xye 



(7) 



Taking the matrix elements of Haso between all three 
t2g orbitals in a similar manner, an effective Hamiltonian 



(2) in these bands, H^go '■^^'^ ^"^ written as 



H 



ASO 
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(8) 



Fig. [2] shows the non-degenerate band structure of 
Hq + Haso for Aaso = 5 meV. From the above Hamil- 
tonian, it may be seen that the lowest energy states and 
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FIG. 2. Dispersion of Hq + Haso with A_b 
Aaso = 5 meV 



50meV and FIG. 3. Band structure corresponding to the Hamiltonian 
Ho + Haso + Ha with Ais = 50 meV, A^ = 10 meV and 
Aaso = 5meV. 



highest energy states mix all three t2g orbitals (with se- 
lected spins), but the middle states only contain dy^ and 
dxz with same spin. 

We now turn to inter-orbital nearest-neighbor hopping 
iJa, which induced primarily by polar lattice distortion 
due to the external electric from inversion asymmetry. 
We can qualitatively understand this as Ti-O-Ti hopping 
process between two neighbor Ti orbitals with different 
parity in z, for example, hopping between d^y — Px ^ dxz 
along y direction. Therefore, the effective for m of H q in 
the basis of ioo orbital bands can be written aa^^I^ 



t2g 



H„ = A, 



X litxi'u 




(9) 



Notice that Ha generates hopping terms from dxy to 
dxz only in the y-direction and from dxy to dyz only in 
the x-direction. Otherwise, the hopping matrix element 
will be an integral of an odd function in a; (or y) and will 
vanish. Fig. [3] shows the spin splitted band structure of 
Hq + Haso + Ha for A^ = 10 meV and /S.aso = 5 meV. 



III. EFFECTIVE SPIN-ORBIT INTERACTION 



Rashba SOI, which takes the form: 



Hr 

AzQ 



I 






-/3iifc_ /32fc_ -^aifc^ 
-/32fc+ /33ifc_ 

-h'^^ -^4fc- 

/34fc+ 

-/33ifc+ /34fc_ /3iiA:+ 

-/34A:- -;3ii/c_ 

(10) 
Here k± = kx i: iky and the order of Ho{k = 0) + Haso 
eigenstates is from highest energy to lowest energy. That 
means (1,0,0,0,0,0) and (0,1,0,0,0,0) correspond to 
the highest energy eigenstates of Ho{k = 0) + Haso- 
Note that there are three energy eigenvalues for i?o(fc — 
0) + Haso, with two Kramers-degenerate eigenstates for 
each. From this Hamiltonian we expect linear in momen- 
tum Rashba SOI (linear Rashba SOI) in bottom and top 
bands and cubic in momentum SOI(cubic Rashba SOI) 
(A ^ ct^k^) in the middle band. The absence of linear 
Rashba SOI in the middle band is due to the fact that 
the middle band at A; ^ contains only dxz and dyz com- 
ponents and, therefore, is odd in z-direction. 



A. Three-Band Model Near fc ~ 

For simplicity, we neglect the fc^ term in the energy 
dispersion in comparison to linear-in-fc terms. In this 
limit, we find the eigenstates of HQ{k = 0) + Haso, and 
then express Ha in that basis. The result is an effective 



The coupling coefficients /S, depend on Aaso and A^;. 
When the band splitting A^; is much larger than Aaso, 
the lowest band is primarily dxy-^ike near fc ^ 0, and we 
can estimate the size of fc— linear Rashba coupling in the 
lowest energy bands ai = A^ a/3i from the second-order 
perturbation, i.e. first-order in Haso and first-order in 



Hn in orbital basis as follows: 



E 



{k,dxy,a'\H'^^^\k,d.j:y,a) = 

{k, dxy, a'\Ha\k', d^^.a") (fc', d^^,(T"\HAso\k, d^y,a) 



^ASO^ 



Ascik) 



Ed^^ik') - Ed^^ik) 

+ [dxz -> rfyz) 

{{(Ty) sm{kxa) - (ax) sm{kya)) 



ai{k X a) ■ z (11) 



where 



AaSqAz 



A, 



ai 



A, 



(12) 
(13) 



This perturbative description of the Rashba SOI breaks 
down at the band crossings of iJo, and we find dramatic 
changes in the strength of SOI as we will see in Sec IV. 
with exact diagonalization analysis. However, we find 
linear Rashba SOI dominates cubic Rashba SOI for small 
k. To see this in more detail, we restore the k^ energy 
dispersion in Eq. ( 10 1 and compute Hr retaining both 
linear and cubic terms. We take A^; = 320 me W^ and 
plot the strength of SOI ai as a function of carrier density 
for several values for Aaso and A2 as we can see in 
Fig. [4] Here ai — Ai^/k^ where A^ is the Rashba 
SOI-induced energy sphtting of the bottom bands. The 
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FIG. 4. ai(meV-A) vs n (lO^^cm-^) at STO surfaces. Up- 
per curve: Az — lOmeV, Aaso = 5meV. Middle curve: 
Az — 5meV, Aaso ~ ISmeV. Bottom Curve: Az = lOmeV, 
Aaso = lOmeV. 

linear Rashba coupling can be identified as the value of 
ai at k — 0. The slope of the plot is proportional to the 
cubic Rashba effect. For various values for Aaso and 
A^, we see that the contribution from the cubic term is 
dominated by the linear term. 




-50 
-100 
-150 

-200 
-250 



FIG. 5. Band structure of 4-band model with Aaso = 9meV 
and Az = 20meV. The first d^y sub-band (dashed lines) is 
assumed to be localized. 



Four-Band Model 



We can extend our theory to four-band model in which 
two dxy sub-bands lies below d^z and dyz band due to 
the strong confinement along 2:-axis(See Fig. Isl). Recent 
DFT calculation^^ shows the first (second) dxy sub-band 
has Ae = 250(50) meV for LAO/STO interface. Most 
of the electrons coming from the polar catastrophe^ 
10^"* cm^^ are localized in the first dxy sub-band as sug- 
gested in density functional calculation of Ref. 129" and 
l30l The main difference with the three-band model is 
that the second dxy sub-band has a much smaller A^ 
since they do not see the large electric field that the first 
sub-band electrons see. Hence, those light electrons do 
not contribute to the anti-localization effect in our pic- 
ture. 



C. Effective Model for Quasi-One-Dimensional 
Channel 

Now we will consider a quasi one-dimensional sys- 
tem at the LAO/STO interface which can be related 
to a nanowire artificially draw n using AFM tip with 
LA0(3.u.c.)/ST0 intcrfac«P12U. We assume there is a 
confinement in y-direction such that wavevector in y- 
direction is quantized as fc" = ^, and degeneracy be- 
tween dxz and dyz bands at /c = is lifted. With quan- 
tized ky, dispersion relations can be written as. 
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(14) 



2mi 



Af 



where AEyiA^z) is the energy splitting due to the y{z)- 
direction confinement. Since the degeneracy is lifted, the 



da;z-like band also has fc— linear spin-orbit coupling, and 
it seems hard to distinguish it from d^^y-like band with 
weak anti-localization measurement. However, we find 
that the relation between chemical potential and spin- 
orbit coupling strength strongly depends on the band. 
As may be seen in FigJ6J if the transport is still domi- 
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FIG. 6. SOI strength versus chemical potential for quasi 
one-dimensional system. The inset shows dispersion of Ho- 
Dashed(Solid) curve is da;y(da;z)-like band. Asy = 30meV, 
Aez = 50meV, Aaso = 9meV and Az — 20meV are used. 
SOI of dxz-like band changes sign at crossing point of dxz and 
dyz band(two solid Hnes in inset). 



nated by da;z-like band, we can see that strength of SOI 
goes to zero at specific value of /io(or fep) because the 
degeneracy between dxz and dyz band is recovered at the 
band crossing point of Hq. When the transport is domi- 
nated by da;y-like band, the evolution of SOI strength as 
increasing chemical potential does not have any nodes. 



IV. EFFECT OF SPIN ORBIT INTERACTION 
ON MAGNETOCONDUCTIVITY 



Quantum corrections to conductivity in Two 
Dimensions 



field as a long-distance cutoff, we derive (relatively) sim- 
ple closed forms of the WAL corrections to the conduc- 
tivity in the limits of (a) purely linear 2D Rashba SOI; 
(b) purely cubic 2D Rashba SOI; and (c) linear quasi-lD 
Rashba SOI. In the first two regimes and in the limit of 
small magnetic fields, we recover the lordanskii, Lyanda- 
Geller, Pikus (ILP) theory ■^'^, which treats the magnetic 
field more precisely by sunrming over Landau levels. In 
the third regime, we obtain similar results to those of 
KettemanrpS. 

We argue that it is essential to use these forms to fit 
WAL data for STO surfaces and LAO/STO interfaces. 
The ILP theory was used by Nakamura et al!^ to deduce 
the cubic /c-dependence of Rashba SOI at the surface of 
STO. (See Eq. (SI) of Ref. 16) However, Cavigha et 
al^ used the Maekawa-Fukuyama (MF) theor}l2l!^ which 
incorporates the SOI simply as a spin relaxation time, 
following Hikami, Larkin, Nagaoka (HLN) theory^ZI. In 
the limit of weak SOI and weak Zeeman splitting, the MF 
theory gives similar results to the ILP theory with dom- 
inant SOI given by cubic Rashba. Therefore, Caviglia et 
al.'s analysis^ could be understood as an indication that 
their results fit the ILP theory with cubic Rashba, as at 
the surface of STO. 

However, since the MF theory and the HLN theory, 
on which it was based, were clearly formulated for a 
physically-distinct situation that is not applicable to the 
LAO/STO interface, it is necessary to compare magne- 
toconductance data to an appropriate theory that takes 
as its starting point either linear or cubic Rashba SOI (in 
2D or ID). We perform such an analysis in the subsec- 
tions that follow and suggest a method for fitting WAL 
at LAO/STO interfaces in the low carrier density region 
which lead us to distinguish the contribution from lin- 
ear and cubic Rashba effects. We discuss experimental 
results from this perspective in the following section. 



The presence of significant SOI changes the universal- 
ity class of the Hamiltonian, and results in a dramatic dif- 
ference in weak field magnetoconductance predictions^Zl. 
Indeed, it is well known that SOI leads to a sign change 
of the quantum correction to conductivity Act. This 
phenomenon, known as WAL, can be used a diagnos- 
tic for the presence of SOI in a conductor. In princi- 
ple, modification of band structure by SOI can be ob- 
served in Shubnikov-de Haas oscillations, but this effect 
will be washed out if impurity scattering is too large. 
The frequency difference between the two Fermi surfaces. 
Aw must satisfy Awr ^ 1. On the other hand, WAL 
is observable under the much less stringent condition, 

2 

f^xx ^ X • I^i'om a detailed fit of the dependence of WAL 
as a function of the density, it is possible to deduce the 
form of SOI, i.e. in our context whether SOI is linear or 
cubic in momentum and, in the former case, whether it 
is due to ID or 2D transport. 

In the subsections that follow, by treating the magnetic 



1. Weak anti-localization due to linear Rashba spin-orbit 
interaction 

Assuming that the dominant contribution to current 
transport is coming from d^y band, we only need to 
consider fc-linear Rashba SOI for small carrier density. 
Therefore, the effective Hamiltonian reads: 



H = h ha ■ fi/fi 

2m; 



(15) 



where a = {ax, cry) is a vector of Pauli spin matrices, 
O/ji = ilfji{sin9, — cos9) is an effective magnetic field 

aik/h and 



of linear Rashba SOI, k = Jkl -\- ky, flm 

ta,n6 — kx/ky. We now consider scattering on short- 
range impurities. For uncorrelated Gaussian disorder 
potential V{x), characterized by the correlation function 

{V{x)V{x')) — 2Tr]/T ^(^ ~ ■^')' quantum corrections to 



the dc conductivity are given by22 
2e2 



Aa 






(27r)2 



CafiPail)- 



(16) 



Here D, v and tq are the diffusion constant, 2D density of 
states and elastic mean- free time, respectively; Cai3pa{q) 
is disorder-averaged Cooperon propagator with a, fi be- 
ing spin indices. Following ILP's approaclP^ESI^ \}^q nia- 
trix equation for the zeroth harmonic of the Cooperon 
propagator which gives the dominant contribution to 
WAL in the diffusive limit, reads 



Here Qmax and gmin are ultra-violet and infra-red cutoffs, 
respectively, with Qmax I/wfTi- If magnetic field is weak, 
one can simplify the calculation by including magnetic 
field as IR cutoff given by q^jj^ — maxjl/Dr^, eS/Zi}. 
(If magnetic field is large, one has to perform the sum- 
mation over the Landau levels^^l). In this paper, we focus 
on weak magnetic field limit, in which case one can ob- 
tain analytical expression for the quantum correction to 
magnetoconductivity: 



Aa(B) - Aa(0) = --—(-A/" + A/^ + A/1 + A/i) 

ZTTh 
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(17) 



It is convenient to rewrite Cooperon propagator in the 
angular momentum basis, in which the singlet J = and 
triplet J — I sectors are decoupled. The eigenvalue for 
the singlet contribution can be readily obtained E'^ = 
D{q^ +q'i) where q'i = 1/Dt^ with r^ being the inelastic 
scattering time. Henceforth, we consider the triplet J — 
1 sector, and find the corresponding eigenvalues. The 
latter requires to diagonalize £ j=i 

Cj^i=Dq'^+— + 2nl^Ti{j^-j^)+ivFTinm{J+q--J- 

— = / W{ip)dip, — = W{ip){l ~ cos{nip))dip, 
To J Tn J 
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+ F+{qso,qB,q4>)- (19) 



Here qg ~ eB/h, and the functions F^{qso,qBTqcji) and 
F+{qso,qB,q<t,) are defined as 
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J^ — Jx ^ ^^y 



q±^ qx± iq 



Here D — VpTi/2 is 2D diffusion constant, W{ip) is scat- 
tering rate for an angle ip, and Ji are vector compo- 
nents of the total angular momentum operator. The SOI 
mixes different components of the J = 1 manifold of the 
Cooperon propagator. By diagonalizing £, one finds 
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where E^ and E}^ are eigenvalues of C corresponding 
to total angular momentum J = and J = 1 sectors, 
and g^Q = 2rL\-^Ti/D = 2a'^ kpTi / h"^ D characterizes the 
strength of SOI. Using these results, one can obtain WAL 
correction to conductivity: 
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(18) 



with qp = 7g^o -t- IGq'i- The magnetic field cutoff q^ 
depends on an arbitrary coefficient ai tha t en ters the 
definition q^ = aieB/h. By comparing Eq. 19 with the 
ILP theory in the limit q^ -^ q% ^ q^^, we find that 
our results coincide if one takes ai — e~^ with 7 being 
Euler's constant. 



2. Weak anti-localization due to cubic Rashha spin-orbit 
interaction 

According to our four-band model, there are two types 
of carriers contributing to the transport. One is the sec- 
ond level dxy band electrons with negligible SOI, which 
we therefore neglect, and the other is middle [dxz + dyz) 
band electrons with fc-cubic Rashba SOI. In this case, the 
effective Hamiltonian for middle band derived from Eq. 



[TOl can be written as 



2m* 



(20) 



where CIr^ — r2fl3(sin36', — cos36') with Hr^ = a^kp/h 
now gives correction to third harmonic of Cooperon, and 
resulting matrix equation for Cooperon in the triplet sec- 



tor can be written as Eq. 17 with 
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(21) 



Now £j=i is readily diagonal in the original basis of J 
and Jz, and its eigenvalues are 
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where q'^^^ = 2n%r^T3/D = 2alk%T3 / h"^ D . Then the 
magnetoconductivity can be written as, 



ACT(B)-Aa(0) = ---(-A/" + A/^ + A/1 + Al]) 
zvr/i 
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= ln 
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(22) 



Again, taking qg=e ''Be/h reproduces HLN theorjEll 
for strong SOI and small magnetic field. 

Clearly, WAL corrections to the magnetoconductance 
are different for linear and cubic Rasliba SOI. Therefore, 
by comparing fits of the experimental data to the above 
expressions, one can try to distinguish between the two 
scenarios. This, in turn, can shed light on the origin of 
superconductivity at the LAO/STO interface. 



B. Quantum corrections to conductivity in quasi 
one-dimensional structures 



In this section, we consider quasi-one-dimensional sys- 
tem confined along the 2/-direction(— ^ < y < ^). In 
this geometry, as previously discussed, the SOI in dxy and 
dxz bands at small carrier density is dominated by lin- 
ear Rashba contribution. We, therefore, concentrate on 
this situation. Given the confinement along j/-direction, 
we need to solve Eq.(17) in real space and impose ap- 



propriate boundary conditions. The singlet component 
of the Cooperon is not affected by SOI, and, thus, the 
corresponding eigenvalue is E^ taken at qy = 0. We now 



concentrate below on J = 1 subspace. The matrix equa- 
tion for the J = 1 components of the Cooperon reads 



^J=iir)Coir,r')^^^6ir,r') 



(23) 
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D 
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Here r = {x, y). The solution of above equation is given 
by. 
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(24) 
(25) 

(26) 



implying zero current in the direction normal to the 
boundary for each spin eigenstates. We note that qy = Q 
does not satisfy the boundary conditions. To find the 
Cooperon propagato r in this case, we first perform a 
gauge transformatioip3I2Zl and simplify boundary condi- 
tions. Let's perform the canonical transformation defined 
by U{y) = e'"'""''^y and introduce \'ipni{r)) = C/(y)|V'm(r)) 
and C{r) = U{y)C,j=i{r)W (y) where C 



C{r) = - + D[{-idxf + i-zdyf] 

-iqso'JxV} 



-2Dq,o[e"'^"'''yjye-'-i--'^y](-idx) 



(27) 



In terms of the new eigenstates, the boundary condition 
reads 



{-idy)\^m{r))\y^±w 







(28) 



and, thus, the zero mode qy = now satisfies the above 
boundary condition. If the width satisfies W <C I/90 
as 1/90 being dephasing length, one can neglect higher 
harmonics, riy > 1, because they are suppressed by a 
factor of Wqrj,. In this regime, the dominant contribution 



comes from 



qy 



— mode. Furthermore, in the 



limit W <^ l/^so, 'C(r) is a slowly varying function of y, 
and can be approximated by its average over y-direction. 
Then, we find that 
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The eigenvalues of L{qx, Q)/D are given by 
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A. Weak Anti-Localization Measurements at STO 
Surfaces 

, Nakamura et ali^ recently reported evidence of a cu- 
bic Rashba SOI at a low carrier density {kpa < 0.3) STO 
surface. They concluded that the bottom dj.j^-like band 
has cubic Rashba SOI. However, as may be seen in Fig. 
H we found that the SOI of bottom dxy-\ike band is dom- 
inated by linear Rashba SOI in the region of small carrier 
density probed in the experiment. 



With the Cooperon propagator in hand, we can now com- 
pute quantum corrections to the conductivity 




At non-zero magnetic field, this expression is modifed by 
introducing an additional cutoff q^ : 



AaioiB) =- 



il + ql 



ql + qB + rqlo 



ql + q% + "^fqlo 

(30) 



c}'3(eV-A^) 
2.5 



A^(meV) 
0.4t 




where r is a width-dependent coefficient that charac- 
terizes the effective strength of the spin-orbit coupling, 
f = (980^^)^/12. The magnetic field cutoff is also 
modified and becomes width-dependent. For weak fields 
B < h/eW'^, the cutoff g| -- e^B'^W'^ /K^ whereas for 
large fields B ^ h/eW"^ , it remains the same as in 2D, 
i.e. q^ ^ eB/h. 

Throughout this section, we have assumed that the sys- 
tem is in the diffusive regime (i.e. all lengths are longer 
than the elastic mean free path) and have derived the 
Cooperon propagator in a quasi-lD system in this limit, 
i.e. assuming that qj^ ^ W ::§> le, where l^ is the mean- 
free path. The results of our calculation can be extended 
to a quasi-lD nanowire whose width is comparable with 
the mean free path lei see, for example, Ref. EH In 
this regime, the magnetic field cutoff should be mod- 
ified due to the flux cancellation effect. In the weak 
{B <C h/eWle) and strong (B > h/eWle) magnetic field 
limit, the cutoff gs becomes g| = 2e'^B'^W^ /{Cifi^le) 
and g^ = 2eBW'^ / {C2hll) , respectively. Here, the coef- 
ficients Ci = 9.5 and C2 = 4.8 are obtained for specular 
boundary condition. 



FIG. 7. Dependence of spin-orbit coupling 03 (eV-A^) 
and spin-orbit energy Ar (meV) on electron density 
n(xl0^^cm~^) in the middle band of STO surfaces. For 
energy spectrum see Fig. [3} Here we used A^ = 35 meV, 
A,4so = 6meV. 

Since this is in contradiction with the data of Naka- 
mura et at^, we assume, instead, that electrons in the 
bottom band are localized (or very poorly conducting) 
and that the observed transport is due to the middle 
band. In that case, WAL should be due to a cubic Rashba 
term since the middle band has only cubic Rashba SOI. 
The calculated values for the cubic Rashba coupling coef- 
ficient as a function of carrier density are plotted in Fig. 

m 



B. Weak Anti-Localization Measurements at 
LAO/STO Interface 

We now discuss SOI using our model for LAO/STO in- 
terfaces and apply it to explain rapidly increasing Rashba 
SOI observed in recent experimentsSiSSl. We consider 
a three-band model, as we did for the surface of STO, 
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and take the values A^ = 50 meV from x-ray absorption 
spectroscopy's and Aaso = 9 meV and A^ = 20 meV 
from the DFT calculations of Ref. [T71 The dependence 
of the strength of SOI, following from the three-band 
model, on chemical potential is shown in Fig. [8J Dashed 
and solid curves correspond, respectively, to bottom (lin- 
ear Rashba SOI) and middle (cubic Rashba SOI) bands. 
Comparison our results with the experiment provides two 
possible explanations for the rapid increase of SOI at a 
specific gate voltage. The first hypothesis is that the 



A^/^,(meV-A) 

A 



/6Q 



-60 -40 -20 




yu(meV) 



FIG. 8. Calculated ai{ABjkx) vs /^(meV). Dashed (solid) 
plots correspond to bottom (middle) band. 

transport at that gate voltage is due to the bottom d^y- 
like band. This band has linear Rashba SOI at A; « 
and a much larger SOI near the avoided crossing with 
the dxz,yz bands. Approaching the crossing causes reduc- 
tion of the gap, which, in turn, leads to an enhancement 
of SOI. For this hypothesis to agree with the observed 
datcP, superconductivity would have to be suppressed 
once electrons start populating the middle band since 
the superconducting transition temperature starts to de- 
crease right after the sharp increase of SOI. It is not clear 
why this should be the case. Furthermore, the observed 
carrier density of the bottom band, as deduced by Hall 
measurements, ^ 10^'^ cm~^ seems to be much smaller 
than the polar catastrophe theory suggests ^ lO^^cm"^. 
Moreover, it is not easy to see why there should be a 
sharp onset density at which superconductivity starts to 
appear. 

The other hypothesis is that these observations are 
dominated by transport from the middle heavy electron 
band. We now apply our four-band model, and assume 
that most of the ^ 10^"* cm~^ electrons predicted by the 
polar catastrophe argument are localized at the interface. 
A much smaller number ~ 10^^ cm~^ of electrons pop- 
ulates the middle band and dominates transport. Elec- 
trons in the second d^y sub-band may contribute to the 
Hall effect, but the WAL phenomena seen in experiments 
would be due to cubic Rashba for small mobile carrier 
density, as at the STO surfaces. This picture also sug- 
gests that superconductivity arises as a result of the ap- 
pearance of electrons in the middle band. 

Yet another possibility to consider is transport in 



quasi-one-dimensional channels, as in wires 'drawn' with 
an AFM tipf^, or appearing spontaneously and are re- 
lated to the formation of tetragonal domains formed in 
the STO below 105 Kpfil. As we discussed above, Rashba 
SOI is always linear in momentum in quasi-lD geome- 
try, regardless of the band index. There is, however, an 
important feature of transport in the quasi-lD case that 
does distinguish between carriers in the different bands. 
At the values of the chemical potential at which the dxz 
and dyz bands become degenerate, the spin-orbit interac- 
tion in these bands vanishes, as shown in Fig. [6j There- 
fore, it would be interesting to revisit WAL data and fit 
the magnetoconductance having a particular scenario in 
mind. 



VI. CONCLUSIONS 

The origin of the physics underlying the ordering 
phenomena of LAO/STO interfaces, namely magnetism 
and superconductivity, is still unclear and controversial. 
However, spin-orbit interaction may provide a window 
into understanding these properties. As we have seen, 
the electrons in the different bands of the LAO/STO in- 
terface have Rashba (i.e. interface-induced) SOI with 
different momentum dependences. Moreover, this depen- 
dence is a strong function of the effective dimensionality 
of the carriers. The momentum-dependence and effective 
dimension are, in turn, reflected in transport measure- 
ments through the dependence of WAL effects on an ap- 
plied magnetic field. In this paper, we have given a simple 
explanation for both the nature of SOI in LAO/STO and 
the surface of STO and also for its WAL signature. 

We find that the sharp increase in the strength of the 
SOI with gate voltage is consistent with conduction that 
is dominated by either the dxy or dxz,yz bands. How- 
ever, the spin-orbit energy scale is predicted to decrease 
at still higher gate voltages. This decrease would steeper 
in the case of the dxy band. Therefore, if it were possible 
to increase the gate voltage until the spin-orbit coupling 
peaks and begins to decrease (as the superconducting 
transition is observed to do) , it would be possible to dis- 
tinguish between these two scenarios. In the quasi-ll? 
case, the difference between the dxy and dxz,yz bands 
may be even more dramatic since the SOI vanishes in 
the latter case at one value of the chemical potential. 
If nominally 2D transport is actually quasi- IZ), as sug- 
gested by recent measurements^Slj then there may be a 
third possible functional form against which WAL data 
on LAO/STO interfaces could be measured. However, 
there are, at present, too many unknowns (such as the 
wire width and spacing) to make a meaningful compari- 
son between theory and experiment. 
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